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Abstract — In this paper we characterize and construct novel 
oversampled filter banks implementing fusion frames. A fusion 
frame is a sequence of orthogonal projection operators whose 
sum can be inverted in a numerically stable way. When properly 
designed, fusion frames can provide redundant encodings of 
signals which are optimally robust against certain types of noise 
and erasures. However, up to this point, few implementable 
constructions of such frames were known; we show how to 
construct them using oversampled filter banks. In this work, we 
first provide polyphase domain characterizations of filter bank 
fusion frames. We then use these characterizations to construct 
filter bank fusion frame versions of discrete wavelet and Gabor 
transforms, emphasizing those specific finite impulse response 
filters whose frequency responses are well-behaved. 

Index Terms — filter banks, frames, tight, fusion, erasures, 
polyphase. 



I. Introduction 

Frames play an increasing role in many signal processing 
applications f\\, ||2l, O. In particular, they have found their 
niche in the engineering community through oversampled filter 
banks (FBs) that implement frame operators |4|, |5|. The work 
presented here is most closely related to |6|, which investigates 
the types of FB frames that are optimal for robust transmission, 
and introduces the concept of strongly uniform frames, proving 
their resilience against erasures. While many constructions of 
oversampled FBs exist, the design of nice oversampled FBs 
— noise and erasure robust FBs with simple reconstruction 
formulas, and built from few-tap filters with good frequency 
responses — remains a nontrivial problem. The aim of this 
paper is to benefit from recent developments in frame theory, 
namely to use fusion frames, to better solve this problem. 
Fusion frames were introduced in Q, where they were called 
frames of subspaces. They have proven useful in different 
settings such as packet encoding [SJ and vector estimation 
from noisy measurements (|9l- In this work, we combine 
fusion frames concepts with FBs, deriving characterizations 
and constructions of FB fusion frames. 

Riding on the success of orthogonal FBs, such as those im- 
plementing discrete wavelet transforms, FB frames emerged as 
a necessary tool to efficiently implement redundant transforms. 
The fundamental works on FB frames H, ||5] characterize 
translation-invariant frames in 1°^ (Z) in terms of the polyphase 
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matrices of the corresponding FBs. In particular, FB frames 
are characterized in [5|, while |4| derives the optimal frame 
bounds of a FB in terms of the singular values of the polyphase 
matrix. One of the most popular examples of FB frames 
is the cosine modulated FB introduced in |10 |. For other 
examples, see |[T], ^ and the references therein. The challenge 
in designing nice FB frames lies in the myriad of constraints 
they must satisfy: to combat noise, their frame bounds should 
be close to each other; to be efficiently implementable, the 
filters should have a low number of taps; in terms of frequency 
response, the filters should have good frequency selectivity as 
well as high stop band attenuation. Even when restated in the 
polyphase domain, the problem of designing such nice FBs 
remains a difficuk subject of interest HIl, El, III21, lfT4ll . 

Frames have become increasingly popular as their redun- 
dancy became a major asset in many applications. This redun- 
dancy is often used to buy stability, resilience to noise and ro- 
bustness against channel erasures. For instance, the optimality 
of frames has been demonstrated in quantization set ups [3 1, as 
well as some source coding schemes. In particular, tight frames 
were proven to be robust to noise, and unit-norm tight frames 
(UNTFs) to be resilient to one or more erasures fTS"], |16|, 
ifTTJI . ifTSl . The generalizations to resilience against packet 
erasures came later with the introduction of a new class of 
frames called fusion frames. Here, individual frame coeffi- 
cients are generalized into projections onto subspaces. Fusion 
frames have already found applications to packet encoding, 
with their robustness against erasures being studied in |8|, 
[9|. Indeed, in ||9l, it is shown that the class of fusion frames 
which, for random signals, are optimally robust against noise 
and erasures correspond to optimal Grassmannian packings, 
termed equidistance tight fusion frames. Fusion frames have 
also been used in other applications, such as modeling feature 
extraction |fT9]| . EOl and sparse recovery fT\\. 

In this paper, we introduce and study a new class of FBs 
named _^/fer bank fusion frames (FBFF). Of particular interest 
is the class of tight FBFFs, which we characterize in terms of 
their polyphase representations and show that they correspond 
to generalized UNTFs. We also present ways of constructing 
FBFFs and tight FBFFs. One example of such a design 
involves taking products of paraunitary polyphase matrices 
with UNTFs. This is exactly the same design process that 
was used to construct strongly uniform tight frames in |6|. 
Indeed, strongly uniform tight frames are a special case of 
tight FBFFs; an alternative generalization of such frames, 
namely totally finite impulse response (FIR) FBs, is considered 
in II22I . Explicit constructions of strongly uniform tight frames 
are presented in ll23l . where the authors design these frames 
in a manner analogous to non-downsampled FBs. In fact, 
the notion of FBFFs can be viewed as a generalization of 
redundant discrete wavelet transforms with non-downsampled 
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Fig. 1. An A'^-channel filter bank with a downsampling rate of M. A 
perfect reconstruction filter bank satisfies $3>* = I. In particular, if the frame 
corresponding to <3? is tight, then there exists A > such that <J>(I>* = Al, 



and so we may take <J> : 



r*. 



FBs of integer redundancy, such as the ones presented in 02411 
and various other works; we pursue this vein, using traditional 
design methods, such wavelet and Gabor transforms, to build 
new examples of FBFFs. In particular, we build discrete 
wavelet and Gabor fusion frames, and in each case, design 
FIR filters that have good frequency responses. 

In the next section, we introduce basic, previously under- 
stood concepts of FBs, frames and fusion frames. Then, in 



Section III we characterize FBFFs in terms of their polyphase 
matrices, and provide several basic constructions of tight 
FBFFs. The remainder of the paper is devoted to constructing 
tight FBFFs whose filters have good frequency responses: 
discrete wavelet FBFFs are considered in Section |IV] while 
Section FVl focuses on discrete Gabor FBFFs. 

II. Background 

In this section, we establish our notation, and discuss 
the previously-known results that we will need in order to 
characterize and construct FBFFs. 

A. Oversampled filter banks and frames 

We consider oversampled FBs Uke those depicted in Fig. [T] 
in which each of the N channels is downsampled by a factor of 
A/, where M < N. So as to be more applicable to real-world 
problems, we restrict ourselves to finite-dimensional spaces of 
discrete periodic signals, namely: 



i{Zp) ■.= {y:Z 



y[p + P] =y[p],p^ z}, 



where P is some positive integer. Equivalently, these are 
signals indexed by Zp := {0, ...,P— 1}, the finite group 
of integers modulo P. In particular, given N filters {'Pn}n=o 
in i{'Zi\ip), the corresponding synthesis FB produces a single 
signal of period MP from N signals j/„ of period P: defined 
as an operator, we have $ : [(.{'Lp)]^ -> ^{'Lmp), 
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(1) 



where ^a/ denotes upsampling by A/, and each ?/„ e ^(Zp) is 
the input of the nth channel of the synthesis FB. The adjoint 



of such an operator is $* : £(Z 
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^*X = {iM{v'n*^)}n=0- 



(2) 



where ^^/is downsampling by M, and i^ajjfc] := ((^„[— fc])* is 
the involution of (^„. A perfect reconstruction FB consists of 
two sets of filters {^Pn\n=o ^^'^ {^n}n=a whose correspond- 
ing synthesis and analysis FBs $ and $* satisfy $$* = I. 

Frame theory provides a mechanism for studying the noise 
robustness of these transforms and, for a given synthesis FB 
$, constructing a corresponding analysis FB $*. Specifically, 
letting I be a countable indexing set, a sequence of vectors 
F — {/i}igi lying in some Hilbert space H is a frame for M 
if there exists frame bounds < A < B < oo such that: 



^ii/iP<Ek/'/')I'^ ^11/11''^/^ 



(3) 
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More generally, {fi}iei is a Bessel sequence if at least the up- 
per bound in (|3]l is satisfied; in this case, its synthesis operator 
F : £2(1) ^ H and corresponding adjoint i^* : H -^ i'^{l) are 
both guaranteed to be linear and bounded: 



Fg 






'«/., (F7)(z) = (/,/,). 



(4) 



The connection between frames and oversampled FBs is the 
fact that every synthesis FB ([T]i and its adjoint (|2| are particular 
types of (frame) synthesis operators and their adjoints Q, 
respectively. Indeed, letting H — £{Zmp), I = Zjv x Zp, 
and identifying £^(1) with [£(Zp)]^, one can show that the 
synthesis operator of $ = {T^^^'~Pn}n=o p=o^ namely: 
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is in fact identical to the synthesis FB ([T]), where T is the 
translation operator: [T^ x)[k] :— x[k — k']. That is, the 
synthesis FB $ is the synthesis operator of the Af -translates 
of the filters {^n\n^o- In this setting, the concept of a perfect 
reconstruction FB corresponds to finding a frame F ~ {/j}iei 
and a dual frame F — {/i}iei so that FF* = I, that is: 

I^Y.^Lh)h, v/eH. 

However, not every frame, nor every perfect reconstruction FB, 
is created equal. Indeed, the whole point of frame theory is 
that certain transforms F will be more stable and will perform 
better than others in the presence of noise, as indicated by 
their condition number. To clarify, the frame condition ([3]) is 
equivalent to having Al < FF* < Bl, where FF* is the 
frame operator. 

FF*^J2f^f^- ^5) 

When H is finite-dimensional, this implies that the least and 
greatest eigenvalues of FF* are the optimal expressions for 
the frame bounds A and B, respectively, making F's condition 
number \J Bj A. As such, for a given application, one typically 
tries to construct frames F whose frame bounds A and B 
are as close as possible, subject to design constraints. Here, 
the best possible situation is when F is a tight frame, that 



is, when A = B, implying FF* = AI. Tight frames are 
especially valuable in the FB setting: in general, any frame F 
has a canonical dual F — {FF*)^^F which satisfies FF* — 
I; however, in the FB setting, it is often the case that the 
canonical dual of FIR filters are no longer FIR; tight FBs are 
an exception to this rule, as (p 
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B. Fusion frames 

The purpose of this paper is to investigate FBs correspond- 
ing to a new class of frames, dubbed frsion frames, which 
are particularly suited to certain applications. Such frames 
arose as a generalization of unit norm frames, that is, frames 
{fi\i& for which ||/i|| — 1 for all i g I. Of particular interest 
are UNTFs — unit norm frames that are also tight — which 
provide optimal robustness against noise and certain types of 
erasures |15|, |16|. Intuitively, the value of UNTFs is that, like 
orthonormal bases, they provide isometric decompositions of 
signals in terms of components of equal weight. However, 
unlike orthonormal bases, UNTFs are nontrivial to construct: 
the easiest examples are the harmonic tight frames obtained 
by truncating a discrete Fourier transform (DFT) matrix; other 
examples are given in |25|, |26|. From an operator-theoretic 
point of view, the significance of unit norm frames is that each 
of the summands fif* of the frame operator (|5]l is a rank-1 
projection operator. Fusion frames generalize this idea so as 
to permit projections of higher rank. 

To be precise, a collection of orthogonal projection opera- 
tors {IlfclfcgK is a. frsion frame if \t& frsion frame operator 
SfceKn/c satisfies Al < J2keK^'' — ^^ ^°'" some < A < 
B < CO. Fig. |2] depicts an example of a tight fusion frame 
in M^ containing four rank-2 projections. These projections 
are onto the four orthogonal complements of the vertices of a 
tetrahedron. Here, X]fe=o ^^^ ~ ^^' where A can be computed 
as the redundancy of the frame as yl = (4 x 2)/3 = 8/3. 
Though superficially a generalization of the concept of a 
unit norm frame, fusion frames are actually a special case 
of them: for each fc G K, letting F^ = {fk,i}ieik ^^ '^^ 
orthonormal basis for the range of 11^, we classically know 
that Ilfc = FkFi. and so: 



Y.u, = j2F,F^ = j2Y.f'^^^fh 



(6) 
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That is, any fusion frame operator may be regarded as the 
(traditional) frame operator (|5]l of a doubly-indexed sequence 
{/fe.JfcsK.jsifc where each subsequence {fk,t}ieu is orthonor- 
mal. Put another way, a unit norm frame is a fusion frame 
if its elements may be partitioned into subsequences, each 
of which being an orthonormal set. From this perspective, 
any unit norm frame which has some mutually orthogonal 
elements, like many traditional wavelet and Gabor frames, is 
an example of a fusion frame. The purpose of this work is to 
better understand and exploit the mechanisms by which these 
orthogonality relationships occur 

Like traditional frames, the value of fusion frames is that 
they permit redundant, stable and noise-robust decompositions 
of signals. But whereas a frame F — {/i}iei breaks a 
signal / into a possibly unorganized collection of frame 




Fig. 2. A tight fusion frame consisting of four orthogonal projection operators 
{Tlk}k—n over R"^. each of which has rank 2. The four corresponding column 
spaces are the two-dimensional orthogonal complements of the vertices of a 
tetrahedron; the unit disk in each of these spaces is depicted. Any given vector 
/ in M^ may be nonorthogonally decomposed in terms of its projections onto 
these four subspaces: / = |(no/ -I- Hi/ + n2/ -|- Ha/). 



coefficients F* f ~ {{f,fi)}iei^ a fusion frame provides 
a more hierarchical decomposition, breaking / into packets 
{Ilfc/lfegK, each of which is itself a collection of coefficients. 



namely 11^^ 



FkF*, where F*f 



{{fjk.t)heu- For 



example, equal rank tight fusion frames, namely sequences of 
equal rank projections whose fusion frame operators (|6]) are 
scalar multiples of the identity, provide isometric redundant 
decompositions of signals in terms of packets of equal weight. 
Such frames are natural generalizations of UNTFs; though 
their existence has only recently been characterized f26|, they 
are known to be optimal with respect to the erasure of entire 
packets H, Q. 

III. Characterizing Filter Bank Fusion Frames 

This article focuses on when $ corresponds to a frsion 
frame. In this section, we emphasize the special case when 
each of the N channels of $$* is an orthogonal projection 
operator; such frames have been shown |8| to be optimally 
robust against the complete loss of data from any given single 
channel. In this section, we characterize FBFFs in terms of 
their polyphase matrices, and then use this characterization to 
provide several basic constructions of tight FBFFs. 



A. Polyphase Matrices of Filter Bank Fusion Frames 

Here, we combine the well-known polyphase characteriza- 
tions of the orthogonality of translates (TT\ and the frame 
bounds of a FB |4| into a new characterization of FBFFs. 
Along the way, we recast these ideas from the infinite- 
dimensional P{Z) setting into the more applicable, finite- 
dimensional £{Zmp) setting. In particular, we define the Afx 1 
polyphase vector ip{z) of some filter ip e liliMp) as: 



viz)-= 






p-i 

<^(")(z):=^<^[m + Mp]z^P, (7) 

p=0 



where, for any m ~ 0, ... M — 1, (^('"^ is termed the rath 
polyphase component of 93. Since ip\m + Mp] is P-periodic 
in p, these entries are technically cyclic polynomials: they 
lie in the ring C[z]/{z^ — 1) of complex polynomials whose 
exponents are integers modulo P. The polyphase space Vm.p 
is the set of all possible such y(z)'s, and is a Hilbert space 
under the inner product: 
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In fact, using properties of the DFT, one can show that 
the mapping (f \ — > vi^)' ^ discrete Zak transform ll28l . 
is unitary with (^, '(/;)^(z„^) = {(p{z),ip{z))p^,^ for any 
(p,ip ^ 1{'Lmp)- Moreover, this fact, along with the easily 
verified relation T'^^^f 1 — > z^'Pip{z), yields that for any 
cc, </9 e £{Zmp), 



{x,T''P^) 



£(Zj 



(x(z),z-VW)i 



p-1 



4ew 



e ^ ),y>(e ^ ))cMe ^ 



p'=0 



= {J-i(x(e^),¥>(e^))cM}W, (8) 

where "•" denotes the variable argument of a given function, 
and T is the nonunitary DFT: 



p-i 

p'=0 



-27rjpp^ 



Taking DFTs of (|8]l yields the fundamental identity: 

(J-(x,r^^V),(z„,))[p] = (x(e^),v(e^))c^.. (9) 

In particular, ^ immediately yields the finite-dimensional 
version of the well-known polyphase characterization of the 
orthogonality of a system of translates: letting x = ipn and 
ip — ipn, we have that {T^^^(pn}p=o is orthonormal if and 
only if 6o[p] = ((p„,T*^P.^„)^(Zj,,p); by taking DFTs, we see 
that this occurs precisely when the evaluation of (p{z) at any 
z = eT^ yields a unit norm vector in C*^. A similar argument 
shows that the Af -translates (pn are orthogonal to those of </?„/ 
precisely when (Pn{z) and (^„/ (z) are orthogonal in C^^ at any 
such z. Moreover, ([HJ provides a means for calculating the 
frame bounds of the synthesis FB $ := {T^^^(pn}n=o, p=o' 
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In particular, letting Bp be the optimal upper frame bound for 

{'^Je^)}^JoMnC^wehave: 



p-i 



p=0 

< (maxPp)||x(2:) 
= (maxPp)|lx|l^(2; 



Ic**^ 



Mp) ' 



(10) 



with equality holding throughout ( |TO| i precisely when x(eT^) 
vanishes for all p such that Bp is not maximal and, at 
the same time, x(e~f^) achieves the upper frame bound of 
{¥'„(eT^)}„ZQ^ for all maximal P^'s. As such, maxPp is 
the optimal upper frame bound for {T^'^^ipn\n=a p=a'^ ^ 
similar derivation provides its optimal lower frame bounds. 
We summarize the preceding discussion as follows: 

Theorem 1. Given filters {p>ri\n=G '" ^(^Afp)-' 

(i) Forafixedn,{T'^^Pifn\n=o is orthonormal if and only 
if iPn[z) has unit norm in (L at any z = e p . 

(ii) The optimal frame bounds for <i> — {T*^''(/3n}„^o p=o 
in £(1,Mp) are: 

A = niinAp, P ~ maxPp, 

where Ap and Bp are the optimal frame bounds for 

W^A^)}^.Zo'inC''whenz = e'-^. 

We emphasize that the novelty of Theorem [T] lies not in 
the results themselves, as (i) is the finite version of a result 
used widely throughout the wavelet literature, and (ii) is the 
finite version of a well-known result on FB frames, namely 
Corollary 5.1 of |4|. Rather, the novelty lies in the realization 
that by combining (i) and (ii), one may characterize FBFFs in 
terms of unit norm frames whose entries are not scalars, but 
rather, polynomials. In particular, if {'■Pn}n=o satisfies (i) for 
all n, then each channel of the FB $$*, obtained by letting 
'fn = 'fin in Fig. [T] is a projection. That is, in this case, $$* 
is the fusion frame operator (|6]l of the projections: 



n„ = 



p-i 
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p=0 



iT''PpnWPiPnr = ^n * (^mIm {v'n * •))• 



The optimal frame bounds themselves are given by (ii). The 
polyphase-based characterizations of Theorem [T] are best visu- 
alized in terms of the corresponding M x N polyphase matrix 
$(z) whose (771, n)th entry is the mi\\ polyphase component 
of the nth filter: 



$„^„(z) = ¥.!") (z), 
and whose N x M adjoint $* (z) has entries: 



(11) 



p-i 



p=0 n=0 



$;,„(z) := [$™,„(z-i)]* = J2{ipn[m + Mp]rzP. 

p=0 

Indeed, (i) means that the columns of <I>(z) have unit norm at 
each z — ct^, while the bounds Ap and Bp in (ii) are the 
least and greatest eigenvalues of $(z)$*(z) at z = e^^. In 
particular, this implies: 



Corollary 2. A synthesis FB <i> generates a tight filter bank 
fusion frame with each channel of <&$* ~ -^I being a 
projection if and only if its polyphase matrix $(2;) has unit 
norm columns and orthogonal rows of constant squared-norm 



M 



at all z 



We note that such polyphase matrices have previously been 
dubbed strongly uniform tight L6J. The remainder of our work 
is focused almost entirely on the construction of tight FBFFs of 
the form characterized in Corollaryl2J that is, polynomial-entry 
UNTFs (PUNTFs). We make this restriction due to the fact that 
such tight, equal-rank fusion frames are known to be optimally 
robust against certain types of noise and erasures |8|, 19|. This 
restriction also permits us to exploit the existing UNTF litera- 
ture. The two exceptions to this rule are in Sections IV andfV] 
where we consider tight, nonequal-rank wavelet FBFFs and 
nontight, equal-rank Gabor FBFFs, respectively. Nevertheless, 
the following theory represents but a fraction of the FBFFs that 
can be constructed using Theorem [T] Moreover, Theorem [T] 
can itself be generalized: by slightly extending an argument 
of II29I , one can in fact unitarily block-diagonalize <1> in terms 
of $(e~f^), writing this MP x NP synthesis operator $ in 
terms of these P matrices of size M x TV, as well as DFTs 
of size Af , P and MP. 



B. Basic Constructions of Filter Bank Fusion Frames 

In this section, we present basic methods for construct- 
ing PUNTFs, namely polyphase matrices $(z) which satisfy 
Corollary l2j thereby yielding tight FBFFs. Most simply, one 
may let $(2:) be a constant-entry UNTF F; the corresponding 
analysis FB $* breaks the signal x in 1{'Lmp) into P 
nonoverlapping blocks of size M, and applies F* to each. 

Example 3. If we let $(z) be the popular Mercedes-Benz 
frame of A^ = 3 elements in M — 2, namely: 



$(z) 



2 -1 

%/3 



-1 

-%/3 



(12) 



then $* is the 3-channel, 2-downsampled analysis FB with 
(^0 = So, (/3i^2 = \{-5o ± y3(5i). Since $(z) is a UNTF 
— unit-norm columns and orthogonal, constant-norm rows — 
for every z, then by Corollary l2] the FB $$* consists of 
three projections of rank ^ which sum to |l, namely those 
projections onto the even translates of i^q. 'Pi and (^2- 

Such PUNTFs are not applicable to most real-world prob- 
lems, as they provide little freedom for designing filters with 
good frequency responses, and also induce blocking artifacts 
in reconstructed data. Nevertheless, the fact that UNTFs cor- 
respond to constant PUNTFs inspires one to see which of the 
very few known methods for constructing UNTFs extend to 
the PUNTF setting. Indeed, most of these methods will extend, 
provided they do not involve division: one may add, subtract 
and multiply the entries of polyphase matrices in the usual 
manner, provided they all lie in the same ring C[z]/(z^ — 1), 
that is, provided the periods of the corresponding filters are 
chosen consistently; division is trickier, and we avoid it; even 
if one alternatively considers FBs over £^(Z), dividing FIR 



polyphase components often leads to unappetizing infinite 
impulse response filters. Specifically, we have: 

Proposition 4. 

(i) // $o(z) and $1 (z) are M x Nq and M x iVi PUNTFs, 
respectively, then so is their M x (Nq + A^i) union 

(ii) // $0 (z) flMc/ $1 (z) are Mq x iVo and Mi x Ni PUNTFs, 
respectively, then so is their MqMi x NqNi tensor 
product $0(2) ® ^i{z). 
(iii) If'^{z) is an MxM unitary polyphase matrix and $(z) 
is an MxN PUNTF then ■^{z)^{z) is also an M x N 
PUNTF 

The proofs of these elementary results are omitted. Note 
that (i) corresponds to simply stacking the channels of two 
FBFFs with identical downsampling rates into a single, larger 
FB. A similar intuitive understanding of (ii) remains elusive: 
unlike wavelet constructions, there does not seem to be a clean 
expression for the frequency responses of the filters defined by 
^q{z)(E)^i{z) in terms of those defined by $o(^) and $i(z). 
Meanwhile, (iii) was previously used to construct PUNTFs 
in [6J in the special case where $(z) was constant. Indeed, 
this approach allows one to make use of extensive wavelet 
literature of paraunitary matrices, and moreover, of the fact 
that all such matrices can be factored into multiplicands of the 
form (I — uu*) + zuu* where u is a unit vector in C*^ |27|: 

Example 5. Let F be the Mercedes-Benz frame ( [T2| i, and let 
'^{z) be the unitary polyphase matrix of the 2-channel, 4-tap, 
2-downsampled Daubechies FB |30|: 

a + hz^^ 
c + dz^^ - 



*W 



- cz 
az 



-1 
1 



(13) 
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where a, d = 2~3(l±\/3), 5,c = 2-5(3t\/3). As such, the 
set of even translates of the corresponding low- and high-pass 
filters ipQ — aSo+c6i+b62+d63 and tpi — dSo—b6i+cS2~a63, 
form an orthonormal basis for €(Z2p). Meanwhile, by (iii) of 
Proposition H the matrix $(z) built as: 

2 + fez-i -2-1(1 -\/3z-i) d-cz-^ 
2 + dz-^ -2-5(^3 + Z-1) -b + az-^ 

is a PUNTF, implying, by Corollary |2] that the corresponding 
3-channel, 4-tap, 2-downsampled FB is a tight FBFF with each 
channel being a projection. That is, the even translates of 
the corresponding filters {(pQ,ifi,Lp2} form a |-tight frame 
for £(Z2p), with the even translates of any particular (p„ 
being orthonormal. The squared magnitude of the frequency 
responses |(^„(w)p = ISfe ¥'n['^]c^^'^"P of these three real- 
valued filters are depicted in Fig. [3] As the first columns of 
$(z) = 4'(z)F and ^(z) are equal, (po is the traditional 
4-tap low -pass Daubechies filter i/jq. Meanwhile, (pi,(p2 = 
I (—'00 i y/Sipi), being combinations of low- and high-pass 
filters, exhibit poor frequency selectivity. As such, this FB $, 
though a mathematically valid example of a FBFF, is not well- 
suited for many signal processing applications. 

Similar problems befall many 3-channel, 2-downsampled 
tight FBFFs: the Mercedes-Benz frame is, up to rotation and 




Fig. 3. The frequency responses of three filters of the 3-channel, 4-tap, 
2-downsampled tight filter bank fusion frame constructed in Example |5] This 
filter bank meets some, but not all, of our criteria for nice filter banks: being 
a tight filter bank fusion frame, it is robust with respect to noise and erasures; 
moreover, it operates quickly, as each filter has only 4-taps; but while i/jq is 
purely low-pass, ipi and ip2 have poor frequency selectivity. 



unit-scalar multiplications, the only 2x3 UNTF ifTSJI . and 
so any 2x3 PUNTF $(z) can be essentially factored as 
^{z) = \E'(z)F, where 'i'{z) is unitary; the three filters of <!> 
are therefore nontrivial linear combinations of the two filters of 
vp, and as such, are likely to have poor frequency responses and 
no good coverage of the frequency domain. One solution to 
this problem is to increase the number of filters, and, perhaps, 
the downsampling rate. In particular, in the next section, we 
construct 4-channel, 2-downsampled FBFFs: as any 2x4 
UNTF consists of two orthonormal bases, we construct a 2 x 4 
PUNTF by applying Proposition Qi) where $o(^) and $i(z) 
are modulated copies of ( pj] ), see Example iTlbelow. Moreover, 
as such FBs correspond to combining two one-level wavelet 
transforms, we are naturally led to investigate how, in a manner 
similar to that of multi-level wavelet transforms, FBFFs can 
be linked together to produce larger ones. 

IV. Discrete Wavelet Fusion Frames 

In this section, we show how by chaining together the tight 
FBFFs of the previous section, we can produce fusion frame 
versions of unitary discrete wavelet transforms (DWTs). Here, 
it is necessary to generalize the notion of a fusion frame to 
that of a weighted fusion frame |7| in which the summands of 
the fusion frame operator (|6]l are accompanied by nonnegative 
weights. Moreover, as we are chiefly concerned with tight 
frames, and may always absorb the tightness constant into 
the weights, we assume the tightness constant is 1. To be 
precise, letting H be a Hilbert space, we say that a sequence 
of operators {Fk : £{Ik) -^ ^}keK induces a Parseval fusion 
frame with nonnegative weights {ck]keE. if Fk-^k =1 for all 
fc e K (or equivalently, that each 11^ = FkF^ is an orthogonal 
projection) and: 



Y, CfcHfc = Y. ''kFkFl = I. 



fcSK 



fceK 



The following result shows how such fusion frames can be 
combined to create new ones: 

Theorem 6. Let {Fk : i{lk) -^ HjfegK induce a Parseval 
fusion frame for H with weights {cfejfegK, o.nd for each 



fc G K, let {Gfc.ij/gLfc induce a Parseval fusion frame for i.{lk) 
with weights {dA;,;}/eLfc- Then, {FkGk,i}k&.,i&uk induces a 
Parseval fusion frame for W with weights {ckdk^i\k&^,i&^k- 

The proof of Theorem [6] is simple, noting: 

{FkGk,i)*{FkGk^i) = Gl [F^FkCkd = GlilGk^i — I, 
for any fc G K and ^ G L^., while: 

/ ^ / ^ Ckdk,i{FkGk.i){FkGk,i)* 



fcsKiGLfc 



E 

feSK 



Cfe^fcl 2_^ akj'^k,l^k,ljJ^k 



iFki > dk.iGk.i.Gt I ]Fi 
Y CkFklF, 



feeK 
= 1. 

The significance of Theorem [6] is that it permits us to build 
intricate FBFFs by composing the operators of simpler ones, 
namely discrete wavelet FBFFs. Though the ideas presented 
here can produce DWTs of any number of levels, we, for the 
sake of clarity, consider only two-level DWTs. To be precise, 
let {fn}n=o ^^ ^^^ filters in an Af-downsampled tight FBFF 
over £{Zm2p) of the type characterized in Corollary [2] In 
the parlance of Theorem l6] we have H = ^(Zjv/2p), with 
each Ik = ^MP and each F,. = $„ = {T^Vnj^lo"' 
corresponding to a single channel of the synthesis FB; viewed 
as operators, we have: 

$„ : e{ZMp) -^ i{I'M^p), *ny = ^n * (Tm 2/)- 

As the M-translates of ipn are orthonormal for any single n 
, we have $* $„ = I. Moreover, as |jl = $$* = J2 "^n*n' 
we have that {^n}n=o induces a Parseval fusion frame with 
constant weights c„ = ^. Next, we choose the Gk/s; the 
manner in which we do so determines the tree structure of 
our fusion frame DWT In particular, to build a fusion frame 
version of wavelet packets — in the second level, each channel 
is passed through a copy of the first-level FB — we choose all 
the Gfc ;'s to be copies of $„'s with dj, ; = "^^ If, on the other 
hand, a traditional DWT tree is desired — the 0th channel is 
passed through a copy of the first-level FB, while the other 
A^— 1 channels are untouched — then the Gq/s are copies of 
$„'s with do,; = ^, while the remaining Gk/s essentially do 
not exist; formally, for fc > 1, we let L^, = {0}, Gk.o = I and 
dk,o = 1. A technical point: one may not choose the second 
level filters to be exactly identical to the first, as they lie in 
different spaces, namely (.{Imp) and ^(Zm2p), respectively; 
rather, they should be chosen as their periodizations; for </?„ G 
^(Zjv/2p), consider: 

M-l 

(fin e ({Zmp), ipn[k] := Y 'Pn[k + MPm], 

m=0 

whose polyphase components satisfy c^,\ ^(z) — if," (z) for 
all z = c~f^; as such, if ^{z) satisfies Corollary l2J then 
^(z) does as well; filters with low numbers of taps are 
essentially equal to their periodizations, and we ignore this 
subtle distinction from this point forward. We conclude this 
section with an example of discrete wavelet FBFFs. 



Example 7. As the 2 x 3 PUNTFs discussed in Example |5] 
yielded filters with poor frequency selectively, we focus here 
on 2 X 4 PUNTFs, that is, 4-channel, 2-downsampled tight 
FBFFs $ in which each channel of $$* is a projection. As 
any UNTF of four elements in two-dimensional space consists 
of two orthonormal bases [151 , we consider the PUNTFs given 
by Proposition pfi), namely ^(z) = [ ^q{z) \ $1(2) ], where 
<i>o(^) ^nd $1(2:) are unitary 2x2 polyphase matrices. That 
is, we consider the FBs obtained by stacking two orthogonal 
DWTs. For example, we can let ^a{z) be the polyphase 
matrix ^(z) of the 4-tap Daubechies DWT, as given in 
Example B] Though we may also take $i(z) to be ^(z), a 
more interesting example is obtained by modulating the ip's, 
shifting their frequencies by |. That is, letting ^o{z) = ^(z) 
and $1(2;) 



[ J 'j']*(-z) yields the 2 x 4 PUNTF: 






cz 



bz 



cz 



"^ jc — idz ^ — j6 + jaz ^ 



The frequency responses of the corresponding filters {(pn}n=a 
are depicted in Fig. |4|a), with a diagram of the corresponding 
synthesis FB $ given in Fig. |4|d). The corresponding tight 
FBFF has redundancy 2, decomposing 4Q-dimensional signals 
in terms of 4 signals of dimension 2Q. Note that such a FBFF 
corresponds to a weighted Parseval frame, of uniform weight 
c„ = i. We further note that in this particular example, 
we additionally have that the zeroth and first channels are 
mutually orthogonal, as are the second and third, implying 
the corresponding FBFF actually breaks such signals into two 
copies of themselves; we do not pursue this fact further here, as 
it does not necessarily generalize to higher-redundancy FBFFs. 
The remainder of Fig. |4] depicts two distinct ways in 
which, using Theorem l6] this FBFF may be iterated so as 
to produce more intricate FBFFs. To be precise, the second 
and third columns of Fig. |4] depict FBFF generalizations of 
traditional DWTs and wavelet packets, respectively. Here, 

since ipn * (Ta/ {fn' * y)) = {fn * (tM fn')) * (tj\/ V), 

we may consider the frequency responses of the equivalent 
filter of any channel: </?„_„' :— ipn * {'\m 'Pn')\ one of the 
true advantages of employing such iterated FBs is that the 
frequency-responses of their equivalent filters are predictable: 
Ifn.n'it)? = \'^nit)\'^\^n'iMt)\'^. In particular, if only the 
low-pass channel is again passed through <1>, a la Fig. We), the 
corresponding 2-level transform is equivalent to a 7-channel 
FB whose frequency responses are depicted in Fig.Qb). Here, 
the four channels corresponding to {(fo.Q,(pQ.i,'PQ,2,'Po.3} 
have a downsampling rate of 4, while the three channels 
corresponding to {(pi,<f2,f3} have a downsampling rate of 
2. By Theorem l6j this FB corresponds to a weighted Parseval 
fusion frame which breaks 4(5-dimensional signals into 7 
components, namely its projections onto four Q-dimensional 
subspaces, each weighted by j, and three 2Q-dimensional 
subspaces, each weighted by i. If, on the other hand, each 
of the four original channels is again passed through $, the 
resulting wavelet packet, depicted in Fig. l4tc),(f), is a full tree 
and breaks 4(5-dimensional signals into sixteen Q-dimensional 
projections, each weighted by j. 

The preceding example highlights the great control in 



frequency design offered by the use of modulation. In the 
following section, we focus exclusively on FBs whose filters 
all arise as the modulates of a single filter. 

V. Discrete Gabor Fusion Frames 

In this section, we characterize and construct Gabor FBFFs. 
That is, we consider the fusion frame properties of collections 
of translates and modulates of a single filter; for y £ £{'Lp), let 
{M.^y)[q\ := e~^^ y[q\ by the modulation of y by a factor of 
p. Here, we assume that the number of channels N is divisible 
by the downsampling rate M; this integer-redundancy case is 
significantly less complicated than the noninteger case 1,31] . In 
particular, letting N = MR, P = QR and (/? G ({^mqr), we 
consider the filters {(fnJnJIf^y where tpn[k] = {M'^"ip)[k] = 

2 TT j /S 71 

e M-R (p[k\. That is, we consider the fusion frame properties of 
the system {T^^PX'^'Vlpi^oT^ ni^~^ of redundancy R. This 
corresponds exactly to an 7\/i?-channel, A/-downsampled FB 
whose filters are regular modulates of a single template (p. 

Letting (p{z) be the polyphase vector (|7]i of ip, one may 
easily show that the mth polyphase component of ipn is 
e MR (p(™)(e K^z). This fact immediately implies that the 
rows of the corresponding polyphase matrix ^{z) are orthog- 
onal: letting n = r + Rk, 



m^wi^))r 



R-1 

X ^e- 

r=0 



AI-1 
fe=0 



■,')r 



^(")(e^F2)((^(™')(e^^z))*; 



evaluating the geometric sum, we see that <l>(z)$*(z) is 
diagonal. Its eigenvalues are thus obtained by letting m = m'; 
Theorem [T] then recovers a finite-dimensional version of a 
previously known result [28] giving frame bounds on integer- 
redundancy Gabor systems in terms of their Zak transforms: 

Theorem 8. Given (p e ({"^mqr), the optimal frame bounds 
o/{r*^P7WQ'V}p=oT^nl5~^ are the extreme values of: 



-27rjr 



R-1 
r=0 



over all m — 0, . . . , M — 1 and all z — e Q« with p G Z. 

In light of Theorem IS] we define the M x R Zak matrix 
<j,Zak(-^-) of if to be the first R columns of its M x MR 
polyphase matrix $(z), sans the external modulation factors: 

In particular. Theorem [8] gives that the synthesis FB $ satisfies 
$$* = Al if and only if the rows of ^^^^(z) are of constant 
norm. Moreover, as shown below, TheoremfTlimplies that each 
channel of $$* is a projection if and only if the columns of 
(J) Zak j- 2^ are of constant norm. Interestingly, when $ is tight 
but each channel is not a projection — the rows of $^^^^(2) 
are of constant norm, but the columns are not — each channel 
is nevertheless a sum of R projections. That is, if (p generates 
a tight A/i?-channel, A/-downsampled Gabor FB, then it is 
necessarily a FBFF, consisting of MR^ rank-Q projections: 




■n —TT 




TT —TT 



(t 2) f2 4 




(d) 



(e) 



(f) 



Fig. 4. The tight filter bank fusion frames constructed in Example|7]using Theorem|6] The first column, from top to bottom, presents the frequency-responses 
and synthesis filter bank diagram, respectively, of the 4-channel, 2-downsampIed tight filter bank fusion frame obtained by combining a traditional discrete 
wavelet transform with one of its modulates. The corresponding fusion frame has redundancy two, breaking signals into four half-dimensional components. 
In the second column, in a manner analogous to traditional discrete wavelet transforms, the low-pass channel is passed through the filter bank a second 
time, producing a new tight filter bank fusion frame, which decomposes signals in terms of seven projections: four of quarter dimension, and three of half 
dimension. In the third column, every channel is again passed through the filter bank; similar to wavelet packets, the coiresponding tight filter bank fusion 
frame breaks signals into sixteen quarter-dimensional components. 



Theorem 9. The Gabor sequence {T^^PM'^"ip}^j^o^^ *i^"^ 
of a unit norm ip e ICZmqr)' 

(i) has each subsequence {T*^''A^'^"<y5}pLo~ being or- 
thonormal if and only if: 



M-l 

E 

m=0 



^(")(z)P = l, 



27rjp 

for all z ^e «« ; 
(ii) is a tight frame if and only if: 

r=0 



M 



(14) 



for all m — 0, . . . , M — 1 and all z = c Q-" , or 
equivalently when, for every fixed m, the R-translates 
of yM(p[m + M»] are orthonormal. 

Moreover, in this case, {T^-^^''+'^'^'>M'^"(p}^^o is nec- 
essarily an orthonormal set for any fixed n and r. 



The proof of this result is given in the appendix. A charac- 
terization of tightness closely related to (ii) is given in |[3T1 . 
Note that in the special case where M = R, then for any ip^^'' 
that satisfies ([14]), namely: 



A/-1 



E I^^^He 



m— 



(0) 



-277j 


-z)P = 


= 1, Vz = 


2^ip 


M 


z e^i-fQ 




2jrjm , 


for all m 




J(e 


M Z) 


= 1, • 



(15) 



,M 



letting (f'-"^^z) = ip 
yields a Lp which satisfies all of the conditions of Theorem l9] 
We conclude this section with an example, using Theoremls] 
to construct a 4-channel, 2-downsampled tight Gabor FBFF: 

Example 10. Let Af = i? = 2 and let Q be arbitrary. Though 
the full 2x4 polyphase matrix is: 



$(z) 






-^(i)(z) 



-j^(i)(-z) 



we see that the rows of this matrix are automatically orthog- 
onal, as guaranteed in general by the proof of Theorem [8] As 



such, we need only consider the smaller 2x2 Zak matrix: 



$ 



Zak 



(z) = 



ip<^"){z) 



o(0)( 



ip^"'{-z) 



In particular. Theorem l9] gives that the discrete Gabor system 
{7~2p_y\/fQn(^| V- 3^^ arising from some unit norm ip in 
£{'Z4q) is tight when: 



|^(")(z)P 
|^(i)(z)P 



|^(")(-z)P 
|^W(-z)P 



(16) 
(17) 



Equivalently, this system is tight if both the even and odd parts 
of ip have norm 2^ s and are each orthogonal to their own even 
translates. In wavelet terminology, we need both the even and 
odd parts of ip to be admissible discrete scaling functions. For 
this to hold, the 4-translates of ip are necessarily orthogonal, 
as are the 4-translates of T^ip — each channel is a sum of 
two projections of rank Q — yielding a FBFF of eight rank-Q 
projections. Note that if we further have: 



|^(")(z)r + |^W(,)|^^l 



Wl 



(18) 



then Theoreml9]gives that the 2-translates of ip are orthonormal 
— each channel is a single projection of rank 2Q — yielding 
a FBFF of four rank-2(5 projections. As noted above, one 
way to ensure ( fTSJ ) is to choose any ip'^^^z) that satisfies ( [T6| ) 



and let tp'^^^{z) — ip^^\—z). Seeking the greatest flexibility 
with regards to the frequency response of (p, we, from this 
point forward, do not make the additional restriction ([TSj, 
settling for the weaker fusion properties guaranteed by ( fTS] ) 
and ( [TTJ i. Indeed, inspired by the max-flat method for con- 
structing discrete Daubechies wavelets ll30l . we seek those T- 
degree polynomials tp'^^^{z) and 1^9^^^ (z) which satisfy ([16]) 
and ( [T7| ) whose corresponding filter ip has the flattest possible 
frequency response at the origin. Specifically, we want the first 
T derivatives of: 



2T-1 



T-1 



T-1 



J2 ^Mz" = J2 ^[2^]^'" + E ^[2p + 1]^'"''' 



k=0 



p=0 



p=Q 



to vanish at z = 1. This equates to a linear relationship 
between the coefficients of (^'^"•'(z) and those of ip^^^z): 



T-1 

E 



(2p+l)! 
(2p + 1 - k) 



T-1 



-ip[2p+l] = - J2 



(2p)! 



=rii 



(2p-fc)! 



fiM, 



for all k — 0, . . . ,T — 1. As such, we can use standard 
nonlinear solvers, such as MATLAB's f solve, to attempt to 
find particular values of the independent variables {(y9[2p]} Zq^ 
whose corresponding ipo{z) and ipi{z) satisfy the Gabor 
tightness conditions ([T6| and ( [T7] i. The numerical evidence 
indicates that these equations have a nonempty, discrete set 
of solutions when T is even. A 20-tap max-flat filter ip 
satisfying ([T6]l and ( [T7] i with T = 10 is depicted in Fig. |5] 



VI. Conclusions 

We have characterized and constructed perfect reconstruc- 
tion FBFFs with desirable properties making them stable and 
robust to noise. We focused on the case where each channel of 




Fig. 5. A 20-tap max flat filter ip whose four modulates and even 

translates form a 4-channel, 2-downsampled tight filter bank fusion frame. 
The frequency response of (p and its three modulates are depicted in (a), with 
the filter itself depicted in (b). The 2-translates of ip are not orthogonal, and 
so, unlike the preceding examples, each channel of the filter bank does not 
correspond to a projection. Nevertheless, in accordance with Theorem lolb), 
the 4-translates of ip ai'e necessarily orthonormal. Indeed, each channel of the 
filter bank corresponds to a sum of two projections, for a total of eight quarter- 
dimensional projections. The fact that this particular filter is even appears to be 
a happy coincidence; most lower-tap analogs do not possess such symmetry. 



the FB is an orthogonal projection. We used polyphase repre- 
sentations to characterize FBFFs and proved that tight FBFFs 
correspond to polyphase matrices with unit-norm columns and 
orthogonal rows of constant norm. Additionally, we presented 
various ways of constructing new FBFFs. In particular, we 
constructed tight FBFFs with wavelet-like structures. Finally 
we studied Gabor fusion frames and designed Gabor FBFF 
filters that have good localization in both time and frequency. 



Appendix A 
Proof of the Theorem[9] 

For (i), note that by Theorem[T| the assumption is equivalent 
to having {T^^^^}p=Q^ being orthonormal. When coupled 



with the following relation: 

{M'^'^ip, T^''PM^"(p) 



-{V,T''P^), (19) 



this fact is immediately 

lQB.~l 



equivalent to having each 
{7-A/p_yv{Qn^|g«-i ^gjjjg orthonormal. 

For (ii), the first characterization of tightness ( [T4] i is an im- 
mediate consequence of Theorem [8] and the fact that the tight 
frame constant A of this UNTF is necessarily its redundancy 



10 



R. For the second characterization, plugging (|7| into ([14]) with 
p = r' + Rq gives that tightness is equivalent to having: 



R-l R.-1 Q-i 



1^ 



'r'=0 



g=0 



Lp[m - 



Mr' + MRq]z-'''-^Ae 



R 

M 



2JJP 

e «« . 



for all m = 0, . . . , Af - 1 and all z = e «« . As the R x R 
DFT satisfies T*T = RI, this, in turn, is equivalent to: 



fl-l Q-l 

E|E 

r=0 g=0 



M¥3[m + M(r + i?g)](z 



i?\-9 



= 1, 



that is, to having the i?-downsampled polyphase vectors of 
each \/Mif[m + AI»] having unit norm at all w = z^ = e Q . 
By Theorem [T] tightness is therefore equivalent to having the 
i?-translates of each y/Mip[m + AI»] being orthonormal. This 
fact then immediately implies that the Afi?-translates of ip 
are orthonormal, which in turn, by ([T9|, implies that each 
subsequence {T''^^''+^'''>M'^"'ip}^^Q is orthonormal. 
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